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Abstract: We present constraints on canonical single- field inflation derived from 
WMAP five year, ACBAR, QUAD, BICEP data combined with the halo power 
spectrum from SDSS LRG7. Models with a non-scale-invariant spectrum and a red 
tilt ns < 1 are now preferred over the Harrison-Zel'dovich model [ng = 1, tensor-to- 
scalar ratio r = 0) at high significance. Assuming no running of the spectral indices, 
we derive constraints on the parameters (ris, r) and compare our results with the 
predictions of simple inflationary models. The marginalised credible intervals read 
ns = 0.962lolo26 ^i^d r < 0.17 (95% confidence level). With respect to previous 
analyses, the portion of the 68% c.l. contours compatible with potentials which are 
concave in the observable region becomes even smaller, but the quadratic potential 
model remains inside the 95% c.l. contours. We demonstrate that these results 
are robust to changes in the datasets considered and in the theoretical assumptions 
made. We then consider a non- vanishing running of the spectral indices by employing 
different methods, non-parametric but approximate, or parametric but exact. With 
our combination of CMB and LSS data, running models are preferred over power- 
law models only by a Ax^ — 5.8, allowing inflationary stages producing a sizable 
negative running — O.OGS^q^qS larger tensor-scalar ratio r < 0.33 at the 95% c.l. 
This requires large values of the third derivative of the inflaton potential within the 
observable range. We derive bounds on this derivative under the assumption that 
the inflaton potential can be approximated as a third order polynomial within the 
observable range. 
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1. Introduction 



Recent advances in centimetre and millimetre wave detector technology are ushering 
in an exciting observational window on the physics of the early universe. The cel- 
ebrated 'inflation' model ^ provides a compelling framework for thinking about 
conditions that might have prevailed in the early universe, yet to date there is no real 
consensus on what is the underlying inflationary model. A well-known implementa- 
tion of the idea relies on scalar field models, for which the inflationary predictions 
can, nonetheless, be worked out in detail. Although it is sometimes claimed that 
inflation 'can predict anything', this view misses the point that, notwithstanding the 
impressive successes of the Wilkinson Microwave Anisotropy Probe (WMAP) [^], 
currently there is a real dearth of hard new facts relating to early universe physics, 
and that the predictions of inflation have yet to be put under any real strain. 

As a consequence there is a need for continual monitoring of the parameters 
of a generic inflation model, as constrained primarily by measurements of Cosmic 
Microwave Background (CMB) anisotropies and surveys of the tracers of Large Scale 
Structure (LSS). Once the values of the parameters of inflation have been determined 
then the implications for the physics of inflation can be worked out. In the long term, 
it is possible that an improved understanding of the physics of inflation could put an 
important constraint on the physics of the early universe, just as an understanding 
of the physics of nucleosynthesis has for a long time been an important constraint 
on physics beyond the standard model of particle physics. 

Some constraints on canonical inflation^ incorporating WMAP five-year data 
(WMAP5) have been published by Komatsu et al. [^, Peiris & Easther |TD[ and 
Kinney et al. [|n[], using different combinations of WMAP, small-scale CMB data 
from ACBAR |12|, Baryon acoustic oscillation (BAO) data from [Q, or SNIa data 
from the Supernovae Legacy Survey [jl4|. Since then, some important new data 



sets have been released, like CMB polarisation measurements from BICEP [jT5| and 
QUaD [0 , the galactic halo correlation function derived from Luminous Red Galax- 
ies in the data release seven of the of the Sloan Digital Sky Survey (LRG7), and the 
luminosity distance of Type la supernovae (SNIa) from the SDSS-II survey [O. 



These new data sets have been a stimulus to our analysis at a time when increas- 
ingly detailed measurements of the CMB from WMAP and Planck |18[, as well as 
a number of sub-orbital experiments such as ACT [T^|, SPT QUIET |^T|, SPI- 
DER PolarBear H and EBEX H, are expected over the course of the next 
couple of years. 

In this work, we will use a combination of the most up-to-date CMB data and of 
the LRG7 halo correlation function for updating constraints on single-field inflation- 



^For recent observational constraints on models with non-canonical kinetic terms, see e.g. 
1^, ^. For constraints on canonical models incorporating WMAP5 data, one should add to the list 
above the update of Ref. Q published in the book |^ . 
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ary models. After discussing our methodology in Section ^ we present a discussion 
of the status of the Harrison-Zel'dovich model in Section ^, which is the simplest 
possible empirical model for primordial perturbations. Then, in Section ^, we show 
our results for the simplest class of inflationary models, those that do not lead to 
any significant running of the tilt in the scalar primordial spectrum. In Section ^, we 
take a more conservative point of view of deriving constraints on the observable part 
of the infiaton potential using different methods (non-parametric but approximate, 
or parametric but exact). Our conclusions are summarized in Section ^ 



2. Methodology 



We perform an analysis of the WMAP |2|, 0, ACBAR [0, BICEP [g and QUaD 



p!6[l CMB anisotropy data, combined with the LRG7 halo correlation function com- 
puted in |2^, Our precise choice of the following subsets of data, hard won by 
all of these experiments, reflects our judgments and decisions that were necessary in 
order to avoid double counting of data. 

Firstly, we prefer to use cosmic variance limited (CVL) data from a satellite 



experiment, where available from WMAP which is CVL up to £ = 530 [^, which can 
be judged to be data of the very highest quality. Note that all suborbital experiments 
derive their absolute calibration from cross-calibration off the WMAP temperature 
anisotropics and ultimately from the CMB dipole, which more broadly speaking 
can be taken as an indication that the suborbital experimental teams agree with 
our assessment of the quality of WMAP data. For this reason we remove a) the 
ACBAR band powers with i < 790 and £ > 1950 to avoid overlap with WMAP 
and to avoid contamination from potential foreground residuals, and b) the QUAD 
TE bandpowers 1-4, and the BICEP TT and TE bandpowers which again overlap 
excessively with WMAP. 

A slightly less appreciated second point is the overlap of QUAD and ACBAR. 



QUAD was originally optimized for making CMB polarization measurements 
and so they chose to observe a small region of sky, and one that overlapped with the 
'Boomerang deep field' in order to improve their absolute calibration inherited from 
WMAP. ACBAR also observed the same field, which they called their 'CMB region 8' 
(see Figure 1 and Table 1 of fl^l)- In terms of gaining independent constraints on the 
CMB temperature spectrum, it is rather unfortunate, but perhaps unavoidable, that 
both ACBAR and QUAD observed the same region of sky. In our final analysis, we 
chose to drop the QUAD temperature data because ACBAR have greater statistics 
on the CMB temperature spectrum (as is evident from their smaller power spectrum 
errors bars) having observed a wider area of sky in total, better optimized for making 
temperature (rather than polarization) measurements. 

Unless otherwise specified, all the results of this paper are derived from this 
combination of data. In some particular cases, we also considered the impact of 
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ancillary data such as the luminosity distance of SNIa from W% and the recent 



determination of the Hubble parameter Hq = 74.2 ± 3.6 km s ^ Mpc ^ by ||30 



We use CosmoMC |^ ^ in order to compute the Bayesian probability distribu- 
tion of model parameters. The pivot scale of the primordial scalar and tensor power 
spectra was set to k^, = 0.017 Mpc^^, as recommended by [^], and we have verified 
that this choice of pivot scale is still close to optimum for the combination of data 
used in our analysis (see Section ^). Apart from the primordial spectrum parame- 
ters (or inflationary parameters) described in the next section, we vary the baryon 
density tUb = fihh"^, the cold dark matter density Uc = flch'^, the Hubble parameter 
Hq = 100/z.kms~^Mpc~^ and the reionisation optical depth r. We assume a fiat uni- 
verse, and so the cosmological constant for each model is given by the combination 
l^A = [1 — i^h + (^c)h~'^]. We also assume a CMB temperature Tcmb = 2.725 K 
and three neutrinos with a negligible mass (excepted in subsection |4.4|) . We set the 
primordial helium fraction to yue = 0.248, a value consistent with the predictions 
of standard Big Bang Nucleosynthesis for a baryon density of Wb — 0.0022.^ In 
order to fit WMAP, ACBAR and QUaD data, we use the lensed CMB and mat- 
ter power spectra and we follow the method implemented in CosmoMC consisting in 
varying a nuisance parameter A^z which accounts for the unknown amplitude of the 
thermal SZ contribution to the small-scale CMB data points assuming the model 



of [§7|- We use a CAMB accuracy setting of at least 1.2. We sample the posterior 
using the Metropolis-Hastings algorithm at a temperature T = 2 (for improved 
exploration of the tails), generating eight parallel chains and imposing a conservative 



Gelman- Rubin convergence criterion |^ of R — 1 < 0.01. 
3. Discussion of HZ model 

The empirical model for the spectrum of primordial perturbations proposed by Har- 
rison |4^, Zel'dovich and Peebles |^ (hereafter the HZ spectrum) has a single 
parameter describing primordial perturbations, namely the primordial spectrum am- 
plitude As- It is interesting to compare its goodness-of-fit with the next-level model 
in terms of number of free parameters, namely a model with a power-law primordial 
spectrum with tilt but no tensor perturbations (this model can be motivated by 
low-scale slow-roll inflation). The general mood is that with the increasing precision 
of CMB anisotropy measurements, the HZ spectrum is on the verge of being ruled 
out in favour of a model with a red tilt. 

However, keeping an open mind about the origin of primordial perturbations, 
the HZ spectrum still holds an allure as a possible indicator of to-be-discovered 



^Under the assumptions of standard BBN, the commonly used value of j/hc = 0.24 corresponds 
to ujh ~ 0.01, which is clearly inconsistent with observations. With the rather modest sensitivity of 
present data to j/hc psi, this choice may not lead to significant bias, but for data from the Planck 



satellite the issue will become relevant [36 
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symmetries of the early universe. As a result we have taken a slightly different tack 
by assuming for a moment that the HZ spectrum is not only a viable but a correct 
description of the data, thereby reducing the dimensionality of the parameter space, 
in order to examine in detail some of the cracks in this model that are beginning 
to appear. Using our basic data set (WMAP5, ACBAR, BICEP, QUaD and LRG7) 
we ran MCMC chains for the HZ model and for the tilted model. Introducing a tilt 
decreases the minimum effective chi square by A(— 21n(£)) = 12.6, which clearly 
implies a tension between the HZ model and the data'^. On top of that, the HZ 
model also leads to different preferred ranges for the cosmological parameters, as can 
be checked from Figure |^ and Table ^ In particular, fixing ns = 1 implies a rather 
high baryon density of cub = 0.0239 ± 0.0007 (at 95% c.l.), which is at odds with 
recent constraints from measurements of primordial element abundances within the 
standard Big Bang Nucleosynthesis scenario. Indeed, a conservative bound based 
only on the deuterium abundance yields uy, = 0.0213 ±0.002 (at 95% c.l.) The 
statistical error shrinks by a factor two when ^He measurements are included ffi 



Taking the deuterium bound on uy, as a prior yields A(— 2 ln(£)) = 17.6 - a difference 
exceeding the Aa level. We conclude that the HZ model is now clearly incompatible 
with CMB and LSS data under the assumption of the standard ACDM cosmology. 

We remark that the HZ model also implies a value of the baryon fraction /b = 
i^b/(i^b + i^c) = 0.173 ± 0.005 (at 68% c.l.) which is about a percent higher than 
the corresponding value in the ns model which impliess /b = 0.162 ± 0.005 (at 
68% C.L). It appears though that, notwithstanding observational uncertainties of 
the total mass of clusters and theoretical uncertainties in the physics of feedback in 
clusters, the slightly higher value of the baryon fraction in the HZ model can easily 
be accommodated by the most recent X-ray observations of clusters (see for instance 



1^). This result merely highlights the current, but rather mild, model dependence 
of the baryon fraction. 

Finally we note from Figure |I] that the preferred value of A^z is zero, which 
results from the fact that both ns and A^z change the ratio of large-scale to small- 
scale power in the Ci. Taken at face value this is arguably another point against the 
HZ model. On the other hand the HZ model does prefer a value of Ho that is slightly 
more consistent with value of Hq = 74.2 ± 3.6 km s^^ Mpc^^ from pO |. 



^While this work was being completed, Refs. [|3| |4| appeared and adressed this issue more 
specificany with different methods. In Ref. [Q, the Bayesian evidence ratio between the HZ and 
power-law models is computed for a data set slightly different from ours, using a prior 0.8 < ng < 1.2 
for the power-law model. The authors conclude that there is strong negative evidence for the HZ 
model. In this paper we do not quote an evidence ratio because this would require a very good 
sampling of the distribu- tion tail for V{ns) in the power law model, four sigma away from the 
maximum, and our runs were not designed for this purpose. Instead, Ref. proposes a method for 
reconstructing the primordial scalar spectrum with no underlying theoretical or inflationary prior. 
The conclusion of this work is still that a scale-invariant spectrum is disfavored, but only weakly 
(consistently, when thoretical priors are removed, one reaches more conservative conclusions). 
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0.92 0.94 0.96 0.98 1 

Figure 1: Posterior probabilities of the parameters of the HZ model (red) and the tilted 
model (black). The dotted orange line denotes constraints on from BBN and the dotted 
maroon line represents the HST prior on Hq. 





HZ tilted 


h 

T 

log [lO^Ms] 


0.0239 ±0.0007 0.0224 ± 0.001 

0.114 ±0.007 0.116 ±0.007 

728+°-°27 n 694+0.032 
nins+0 036 n n8'?+0 032 

U.iU»_o.o34 U.U»d_o.029 

3.14 ±0.07 3.12 ±0.06 


0^8 


0.853 ± 0.048 0.822^°:°^^ 



Table 1: Mean parameter values and bounds of the central 95%-credible intervals for the 
HZ model and the tilted model. 



To summarise this section, we have argued that the HZ model faces pressure both 
on grounds of goodness of fit, but also on grounds of poor astrophysical consistency. 
As mentioned later in section [4.4| , the inclusion of neutrino masses would not alleviate 
the pressure on the HZ model. 
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4. Constraints on the simplest models of inflation 



4.1 Motivations 

Most inflationary model discussed in the literature lead to a power-law scalar spec- 
trum, i.e., to a negligible running of the scalar tilt. This follows from a well-known 
argument which can be summarised in the following way. The function accounting 
for the Hubble parameter as a function of the field (or as a function of the number 
of e-folds) can be expanded around any value in an infinite hierarchy of slow-roll 
parameters e„, each one accounting for the logarithmic derivative of the previous 
term e„_i (see for instance ^ or Appendix |A.2|) . Current bounds on the 



tensor-to-scalar ratio r and scalar tilt ns imply that the first two parameters ei and 
€2 are much smaller than one, as expected during slow-roll inflation. So, a sizeable 
running as on observable scales can only be generated when one of the parameters 
e„ with n > 3 is roughly of order one or larger when observable scales leave the 
Hubble radius during inflation, i.e., when some derivatives of the first two slow-roll 
parameters are very large. If this is not the case, the running will be of the order of 
€162, i.e., at most of the order of as ~ 0(10""^). 

On the other hand, the number of e-folds required between that time and the end 
of inflation (which depends on the scale of inflation, on the details of inflation ending, 
and on the efficiency of the reheating mechanism) can be conservatively assumed to 
lie in the range from 30 to 60 e-folds. For any smooth potential, a sizeable running 
|as| > 0.01 implies such large derivatives of the first two slow-roll parameters in the 
observable range that infiation would end very few e-folds after galaxy scales leave 
the Hubble radius, not even reaching 30 inflationary e-folds [0, . This argument 
can only be evaded: (i) with very special potentials (incoporating sharp features, or 
such that their Taylor-expansion involves very high-order coefficients, see for instance 
[0, ^ H^); (ii) within set-ups involving several inflaton fields, leading to a phase 



transition or even to several short stages of infiation, such that the required 30 to 60 
infiationary e-folds are not contiguous. Since these situations are beyond minimal 
requirements, and since the sensitivity of current data to running (Aas ~ C(10~^)) 
is far too small for probing generic slow-roll model predictions (as < (9(10~^)), it 
appears sensible to perform an analysis where the running is assumed to be negligible, 
in order to preserve the simplicity of the infiationary paradigm. We will relax this 
theoretical prejudice in the next section. 

We thus consider all infiationary models which can be described by the pri- 
mordial perturbation parameters consisting of the scalar amplitude and spectral 
index {As,ns), and the tensor-to-scalar ratio r (both defined at the pivot scale 
fc* = 0.017 Mpc~^). In single-field infiation, deep in the slow-roll limit, the ten- 
sor spectrum shape is not independent of the scalar one. We will consider a tensor 
spectrum with a tilt = —r/8, as predicted for canonical single-field infiation at 
first-order in slow- roll. 
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4.2 Results for the basic data set, and implications for inflation 

In Table we present constraints on each parameter of this model, using only CMB 
data in the first column, and our reference data set (CMB plus LRG7) in the second 
one. In Figure |^ we focus on the joint probability in the (ns, r) plane, for the reference 
data set. Compared to the results obtained by @, our constraints are shifted towards 
slightly lower tilts and significantly lower tensor-to-scalar ratios. The marginalised 
95%-credible interval for the tilt is given by = 0.9621qq26 (to be compared with 
the result of the previous section, no, = 0.9551q;q26' obtained for low energy scale 
inflationary models with r ^ 0). At the same confidence level, our result for the 
tensor-to-scalar ratio is r < 0.17. This corresponds to an upper bound on the energy 
scale of inflation 



3^3 r rrip 
128 



< 8.9 X 10 



-12 



rriT 



(2.1 X 10'*'GeV) 



(4.1) 



at 95% C.I., where mp = G^^^^ stands for the Planck mass (the reduced Planck mass, 
referred to in the next sections, will be denoted as Mp = (SnG)^^^'^). Bounds on 
quantities which are not sharply constrained by the data, like r or Vq, depend mildly 



on the choice of parametrisation, as illustrated in [55|, but for reasonable choices this 
will not affect our qualitative conclusions. 





CMB only 


CMB±LRG (DR7) 


CMB±BAO (DR7) 


h 

T 

log [lO^O/ls] 
r 


0230+°-°°^3 
n 746+0-066 

n nqi 

U.Uyi_o.o34 

3.08 ±0.08 
n q7ci+o.039 

< 0.29 


0.0227 ±0.0011 

Q 1 1 C-+0.008 

0.083l°;°i 

q 1 9+0.06 
"-"•-■-^-O.OT 

9fi2+°-°2^ 
< 0.17 


0.0226 ±0.0011 
0.113 ±0.007 
0.708 ±0.030 

U.UOd_o 030 
q 1 1 +0.07 
'-'•-'-^-0.06 

961+°-°27 

U.»Ul_o 026 
< 0.18 



Table 2: Mean parameter values and bounds of the central 95%-credible intervals for 
the parameters of the vanilla-l-r model and various combinations of data sets. For the 
tensor-to-scalar ratio r the 95%-credible upper bound is quoted. 

As pointed out in previous works (see for instance [^, 0), current data mainly 
raise some tension for inflationary models with a concave potential in the observable 
range. Most of the 68% allowed region lies below the convex potential limit, and 
convex models with a red tilt in the range [0.93-0.99] are comfortably allowed by the 
data. These include for instance: 



'Natural inflation' 
the prediction ns 
(95% C.I.). 



56 



with the potential V{(j)) = A'' [1 ± cos(0//)] leading to 
— mp/{8TTp), which is allowed only for / > 0.73 mp 
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Figure 2: Marginalised joint 68%-, 95%- and 99.7%-credible contours in the (ns, r) plane 
for our reference data set (CMB plus LRG7). The orange line shows the limit between 
inflationary models with a concave/convex potential in the observable range. We show in 
black the regions corresponding to popular concave models with a quadratic, quartic or 
exponential potential. 



More generally, 'new inflation' type models where the inflaton rolls away from 
unstable equilibrium with a slope given by V{(f)) = Vo[l — For a = 2 

and assuming |0| <^ /i when observable scales leave the Hubble radius, our 
constraints imply 1.5 < fi/rrip < 4.0 (95% c.L). Inflationary models with 
a > 3 are all in good agreement with observations. 



The symmetry-breaking potential V{(f)) = A(0^ — 0q)^ |]57| , |58| is convex in its 
observable region for 0o < 6.6 mp. We find that this model is in agreement 
with observations for 0o > 2.1 mp and A ~ (9(10~^^). 

'Hybrid inflation' with a logarithmic slope (caused for instance by one-loop 
corrections during global SUSY inflation ) , which predicts a small r and ns 
close to 0.98. 

Inflationary models predicting a red tilt for scalar perturbations and a tensor- 
to-scalar ratio ~ 0{{ns - if), as R + R^/{QM^) [g or Standard Model Higgs 
inflation [^, are in good agreement with observations (as well as models with 



an even lower r, see e.g. Refs. |^ p2|). 
Among concave models, we can single out a few generic models for further discussion: 
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• Inflation with a quadratic potential V = ^w?(j? (often dubbed 'chaotic infla- 
tion' after Ref. [^) lives along the line r ~ — 4(ns — 1) during slow-roll. The 
position of the model along this line depends on the number of e-folds A^* be- 
tween the time of Hubble crossing for the pivot scale, and the end of inflation. 
In turn, this number depends on the reheating temperature. Due to uncertain- 
ties on the reheating stage, we consider here a plausible range 50 < A^* < 60, 
for which quadratic inflation is still within the 95%-credible region in the (ns, 
r) plane, as can be seen in Figure |. Note that updated constraints on the 
mass can be inferred from the combination {?)'kA^Y^'^/{2N^) which is equal to 
(m/mp) in this model. Varying As in the 95%-credible range given by Table 
0, and A^* between 50 and 60, we find that the mass of chaotic inflation models 
lies in the range from 1.2 x 10~^mp to 1.5 x 10~^mp. 



Inflation with a quartic potential V = A0^/4 lives along the line r = — 16/3(?t,s — 
1) during slow- roll. For this model, there are indications that the number of 
e-folds should be increased by four, for the reasons summarised in ||6^. In 
Figure ^ we show the line corresponding to quartic inflation for A^* < 64. This 
model is found to be outside of the 99.7%-credible region in the (ns, r) plane. 



• Inflation with an exponential potential V = exp [-^/p/2{(p/Mp)] is called 
power-law inflation [Q, because the exact solution for the scale factor is given 
by a(t) oc t^. This model is incomplete, since inflation would not end without 
an additional mechanism which stops it. Assuming that such a mechanism 
exists and leaves unmodifled its predictions on cosmological perturbations, we 
can constrain this model since it predicts r = — 8 {ns — 1 ) and ns = l — 2/{p — l). 
Only models with 78 < p < 246 are found to lie within our 95%-credible region 
in the (ns, r) plane. 



Hybrid or 'false- vacuum' inflation with a quadratic slope and an effective po- 
tential V = Vo + ^m'^cj)'^ is also compatible with the data provided that Vq is not 
too large with respect to the quadratic term. In this model, the two slow-roll 
parameters are related through ai]"^ + e — rj = 0, where a = 87rVo/(m^mp). 
For a > 80 this relation is not compatible with our 95%-credible region in 
(ns,r) space. We conclude that independently of the mechanism responsible 
for the end of hybrid inflation, this class of models should obey to the con- 
straint Vq < Sm'^mp. Note however that this model is usually invoked as a way 
to obtain ^ Mp, and to avoid large radiative corrections during inflation 
[3^1. This can be achieved when the term Vq dominates over the quadratic one, 
i.e., in the large a limit - a situation constrained by current data. 
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4.3 Impact of extra data sets 



4.3.1 Baryon acoustic oscillations 

The results presented in the third column of Table ^ show that for the model consid- 
ered in this section, the constraints on the primordial spectrum parameters do not 
degrade significantly if we replace the LRG7 power spectrum by constraints on the 
baryon acoustic oscillation scale derived from the same survey data This is a 



clear indication that the shape information contained in the galaxy power spectrum 
data is not very relevant here - the improvement over the results from a CMB only 
analysis stems entirely from the geometrical information; it also implies that our 
estimates are unlikely to be affected by uncertainties in the modeling of non-linear 
structure growth at small scales. 

4.3.2 Supernova luminosity distances 

Even though the luminosity distances of type la supernovae are not a direct probe 
of the primordial perturbations, their ability to constrain the expansion history of 
the low-redshift universe may still contribute useful information about inflation by 
helping alleviate parameter degeneracies in the model. The most recent compilation 
of SN data was presented by Kessler et al. , who give two different sets of lumi- 
nosity distances, derived from the same observations, but using either the mlcs2k2 
or SALT-ll lightcurve-fltting algorithms. In order to test whether SN add any useful 
information to our analysis of inflationary parameters we have separately combined 
the two luminosity distance data sets with the CMB+LRG7 data. In Figure |^ we 
plot the resulting constraints in the (ns, r) plane. The discrepancy between the 
results for the two light curve fitters is evident, and has already been pointed out 



in Ref. JT^ in the context of the dark energy equation of state parameter. Whereas 
adding the salt-ii data has no appreciable effect on our and r constraints, the 
mlcs2k2 data would lead to a significantly tighter bound on the tensor-to-scalar 
ratio. This effect can be traced back to ]V[LCS2k2's apparent preference for large 
values of Qm- Given that at the moment the SN data fail a basic self-consistency 
test between different methods, we choose to refrain from combining them with the 
more robust CMB+LRG7 data in the following sections. 

4.3.3 Direct constraints on the Hubble parameter 

We studied the impact of imposing the recent determination of Hq from the measure- 
ment of nearby supernovae with the Hubble Space Telescope , which corresponds 



closely to Hq = 74.2 ± 3.6 km s~^ Mpc~^ for the ACDM cosmology (we verified that 
'post-processing' the Markov chains with a Gaussian prior on Hq yields virtually 
identical results to those obtained using the weakly cosmology dependent Hq likeli- 
hood code provided by [0| and included in the October 2009 version of CosmoMC). 
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We found that since CMB+LRG7 prefer a slightly lower value of Hq at la, then 
imposing the Hq constraint shifts ub and hence ns to slightly larger values. 



4.4 Impact of neutrino masses 

In this subsection, we address the probably most well-motivated extension of the 
basic cosmological model, massive neutrinos. In pB[, the inclusion of a neutrino 



mass parameter was shown to weaken parameter constraints on ng and r enough 
to affect conclusions about the A(/)^ model of inflation in particular, mainly due to 
a parameter degeneracy with r. Interestingly we find that the more recent data 
effectively break this degeneracy: the constraints for the vanilla+r-|-m,y model in the 
(ng, r) plane are essentially identical to the ones of the massless neutrino model, see 
top right panel of Figure ^ We find an upper limit on the sum of neutrino masses 
of '^rriu < 0.64 eV (at 95% c.l.), consistent with the results of a recent analysis 
with similar data sets One might also wonder whether the presence of massive 
neutrinos could help alleviate the problems of the HZ model by virtue of the free- 
streaming-induced suppression of the matter power spectrum at small scales. This 
does not happen to be the case here, however; neither the baryon density nor the 
spectral index are significantly shifted with respect to the massless neutrino model. 

4.5 Impact of recombination uncertainties 

One of the remaining major theoretical uncertainties in the calculation of CMB 
angular power spectra is the physics of recombination. The standard approach relies 
on an effective description of hydrogen recombination assuming a 3-level hydrogen 
atom, as implemented in the widely-adopted RECFAST code [|7D|, |7T|, whose latest 
version 1 .4.2 J7|, is the default model used in CAMB. The authors of RECFAST 
estimate that the computed ionisation fraction Xc{z) is accurate at the percent level. 
However, any error in Xc{z) will inevitably propagate to the Ces, affecting in particular 
the small scales, and can potentially lead to biased parameter estimates. 

Recently, there has been a lot of progress in understanding the physics of recom- 



bination (see for instance Refs. |7^, |75|), and it seems that the accuracy required for 
Planck-quality data has now been reached. The results of the latest state-of-the-art 
recombination calculations can be implemented in a parameter estimation analysis 
with the help of RICO [0, an interpolation code for Xe{z). In order to estimate the 
bias from recombination errors on the inflation-related parameters, we performed 
an analysis with RICO (including a 'fudge function' accounting for corrections due 
to Raman scattering and Ly-a radiative transfer effects, see Figure 1 of Ref. fT^ ) 
instead of RECFAST. 

Our results are shown in the bottom panel of Figure ^; the credible regions in the 
(rig, r) plane are shifted slightly towards bluer tilts, though not at an alarming level. 



We confirm the conclusion of Ref. ||7^ that for present data, the current version of 
RECFAST is sufficiently accurate. 
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Figure 3: Marginalised joint 68%- and 95%-credible contours in the (ns, r) plane, il- 
lustrating the impact of various systematic effects. The dotted black lines represent the 
results of our 'default' setup in all three plots and are identical to the contours presented 
in Figure |2|. Top left: Results for CMB-I-LRG7 data alone, and combined with luminosity 
distances from the supernova compilation of Ref. analysed with the mlcs2k2 (solid 
purple lines), and salt-ii (solid pink lines) lightcurve fitters. Top right: Comparison of 
results for the massless neutrino and massive neutrino (orange lines) cases. Lower left: 
CAMB's default recombination code RECFAST 1.4.2 compared to recombination calculated 
with RICO, including a 'fudge function' which accounts for corrections due to Raman scat- 
tering and Ly-a radiative transfer effects (red lines). Lower right: Results for CMB-I-LRG7 
compared to CMB+LRG7 combined with the Hq measurement from |3^] (solid maroon 
lines). 

5. Conservative constraints on the observable inflaton potential 



5.1 Motivations 

The results of the previous section assume a negligible running of the scalar and 
tensor tilts. Indeed, these runnings are completely generic in inflation, but typically 
as small as as < 0(10^^) for the simplest slow-roll models, as already stated in 
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Section ^ In other words, we implicitly supposed in Section ^ that the inflaton 
rolls down a smooth potential well satisfying the slow-roll conditions and providing 
enough inflationary e-folds for matching on to the post-inflationary expansion era. 
These assumptions maximise the beauty and simplicity of the inflationary paradigm. 
However, nature sometimes prefers some degree of complexity, and in principle infla- 
tion could take place along a complicated potential, or in a multi-dimensional field 
space, or in several stages associated to distinct mechanisms. Even in that case, 
inflation remains a key ingredient in the cosmological evolution, necessary to solve 
the flatness and horizon problems and to generate primordial perturbations. If we 
take the point of view of studying what the data really tells us on the inflaton po- 
tential (with the only assumption that when cosmological scales leave the Hubble 
scale, inflation is driven by a single scalar field with a canonical kinetic term), the 
best we can do is to measure the primordial spectrum or the underlying inflaton 
potential within the region probed by cosmological observations, with no attempt to 
extrapolate beyond that range. 

It is sometimes argued that fitting the data with a large running in the scalar 
spectrum (of the order of as ~ O{10^^) or more) is not consistent with the infla- 
tionary paradigm, because - as already argued in section ^?T] - such a large running 
signals that inflation ends soon after observable scales leave the Hubble radius. It 
is worth clarifying this argument. Of course, as long as the purpose is to constrain 
inflation self-consistently, one should check that any large running model providing 
a good fit to the data is indeed compatible with the assumption that the expansion 
is accelerated as long as observable scales leave the Hubble radius. If this is the 
case, such a model cannot be eliminated without invoking some criteria of simplicity. 
Indeed, the observable part of the potential can always be extrapolated in such way 
to accommodate an arbitrary number of e-folds. For models with a large tilt run- 
ning in the scalar spectrum, the simplest extrapolation schemes would suggest that 
inflation ends soon after cosmological scales leave the Hubble scale, but it is always 
possible to design the potential (i.e., to introduce enough derivatives) in order to 
extend the duration of inflation as desired. Alternatively, and as already mentioned 



in section it is always possible to stick to the assumption that the potential is 
smooth and inflation ends very quickly, provided that later on, one or several extra 
inflationary stages (associated with other inflatons) sum up to the desired 30 to 60 
e-folds of accelerated expansion. This situation can even be argued to be generic in 
some particle physics frameworks see e.g. [0). 

5.2 Methods 

The shape of the inflationary potential within the observable region can be con- 
strained through various approaches that make accurate inflationary predictions: 

• The scalar primordial spectrum can be Taylor-expanded at various orders (in- 
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eluding a running of the tilt, a running of the running, etc.). Then, the free 
parameters of the model consist of the scalar spectrum amplitude, its various 
logarithmic derivatives at the pivot scale, and the tensor-to-scalar ratio. For 
consistency, the tensor spectrum shape should be fixed by the hierarchy of 
self-consistency conditions truncated at some order (see Appendix |A.1| ). The 
data are used to provide constraints on these parameters {As, ns, as, r). 
A posteriori, these constraints can be converted into constraints on the scalar 
potential derivatives (or on combinations of these derivatives called the 'po- 
tential slow-roll parameters') by making use of some analytical formula valid 
at a given order in slow-roll (with the caveat that slow-roll can be marginally 
satisfied in this context). This approach was chosen for instance by 0. 

One can choose a parametrisation of V{(j)) in the observable range (if the 'ob- 
servable potential' is assumed to be smooth, a Taylor expansion is adequate). 
For each potential parameter set, the primordial spectra can be computed an- 
alytically or numerically and the model can be fitted to the data. In this case, 
the data provide direct constraints on the potential parameters, without rely- 
ing on any slow-roll expansion in the case of a numerical computation. The 
result depends on the initial value of 0. However, if one does not assume that 
inflation started just before the observable region is crossed, there is a unique 
choice for 0, corresponding to the inflationary attractor solution in phase space. 
This approach can be easily followed by making use of the CosmoMC inflation 
module^ released together with This module performs a fully numerical. 



accurate and fast computation of the primordial spectra, based on the integra- 
tion of background and perturbation equations during inflation (see Appendix 



E3 



In order to relax the assumption on (p, one can perform a similar analysis tar- 
geting the function H{(j)) instead of V{(f)). Since each H{(j)) corresponds to 
a unique V"(0), this method is equally appropriate for constraining V{(j)), and 
naturally incorporates models for which inflation starts just before entering the 



observable region. This approach was followed in I^O], ^ ^ ^ ^ |10|, with 
various analytical or semi-analytical schemes for the calculation of the primor- 
dial spectrum (based on an expansion of H{(f)) into Hubble Slow-Roll (HSR) 



parameters). It was also followed in |^2|, |8^, |^ with a numerical calculation 



of the spectrum for Taylor-expanded functions H{(j)). Although reconstructing 
H[(j)) is slightly more general than V{(f)), we will not follow this approach in 
the current paper because constraints obtained directly on the potential pa- 
rameters (for instance on V /V instead of H'/H"^) are slightly more convenient 
and suggestive for inflationary model-building. 



*http : //wwwlapp . in2p3 . f r/~valkenbu/inf lationH/ 
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One can try to reconstruct the function H{N) where = Ina is the number 
of e-folds, using its expansion in Horizon Flow Functions (HFFs) evaluated at 
the time when the pivot scale crosses the Hubble radius during inflation. 
Combinations of the first three HFFs at the pivot scale can be related to the 



primordial spectra using the second-order slow- roll approximation The 
same approximation also provides relations with the derivatives of the infiaton 
potential (see Appendix |A.2| ). This approach was chosen in many papers, for 
instance ||8^, |86[| . 

In this work, we will mainly compare the first and second approaches, although 
including also the fourth one for comparison. Hence, we will first fit the data with 
primordial spectra described by the set of parameters {As, ns, as, r), and use second- 
order slow-roll formulae to translate our results either in terms of convenient combina- 
tion of the potential parameters {V, V, V", V"), or in terms of the HFF parameters. 
Next, we will check whether these results agree with a direct reconstruction of the 
HFF parameters (using the second-order slow-roll apprimation to relate HFF param- 
eters with the primordial spectra). Finally, we will relax any slow- roll assumption 



and constrain directly the same quantities using the numerical module of |78], as- 
suming that the observable potential V{(p) can be described by an order two or three 
Taylor-expansion. 

5.3 Results 

5.3.1 Flat priors on the spectral parameters including running 

We first run CosmoMC with flat priors on the spectral parameters (In As, ns, as, f')- 
For each parameter, the mean values and marginalised bounds on the central 95%- 
credible interval are given in Table ^ for three cases: with CMB data only, with 
our default dataset CMB+LRG7, and for the same dataset when r is kept fixed 
to zero (low-scale inflation limit). As in previous papers on the subject (see for 
instance [^, ^), a large negative running is preferred, since the value as = is 
always above the 95%-credible interval. Adding the LRG7 data leads to a slightly 
smaller mean value as = —0.063, but does not increase the level of significance at 
which running spectra are preferred. With the CMB+LRG7 data set, introducing a 
running decreases the minimum effective chi square by A(— 21n(£)) = 5.8, showing 
that a non-zero running is preferred by the data, but not with a high degree of 
significance^ (we recall that between the HZ and power-law model the effective chi 
square decreases by twice the same amount). In Figure ^ (red curves), we show the 
marginalized likelihood contours for two-dimensional projections of the parameter 



^ While this work was being completed, Ref. |43 appeared, based on a slightly different data 



set. By computing the 'Bayesian evidence' ratio between the running and power-law models (with 
a top- hat prior on as in the range [-0.1,0.1]), this analysis concludes that 'running is not disfavored 
by the data nor required in modeling the data'. 
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Figure 4: Constraints from the CMB+LRG7 data set on combinations of the HFFs 
evaluated at the pivot scale (ei, 62, 6263) assuming: (red) €3 = and the first-order slow- 
roll approximation for the computation of the primordial spectra; (black) €3 7^ and the 
second-order slow-roll approximation. 



space (ng, as, r). As usual, the running is found to be slightly correlated with the 
tensor-to-scalar ratio, and for large values of r the data are compatible with more 
running. However, in the small-scale inflation limit (i.e., when we fix r = 0), running 
is preferred with roughly the same level of significance as in the r 7^ model, since 
in that case we find as = -0.046lo:o39 at the 95% level. 

In the running model, the constraint on r degrades by a factor 2 with respect 
to the power-law model. The allowed range for rig enlarges by a factor 1.5, but the 
mean value remains identical. Indeed, we choose our pivot scale at k^, = 0.017 Mpc^^ 
in order to remove most of the degeneracy between and ag, as can be checked in 
Figure ||. 

5.3.2 Flat priors on combinations of the HFFs at the pivot scale 

Next, for comparison, we run CosmoMC with fiat priors on combinations of the first 



three HFFs evaluated at the pivot scale, namely ei, 62 and €263 (see Appendix |A.2 
for definitions). We use the second-order slow-roll approximation in order to relate 
these parameters to the spectral parameters (ng, ag, r). In Figure ^, we show our 
marginalized likelihood contours in two-dimensional projections of the parameter 
space (ei, €2, e2e3)- In Figure |, we compare our results with those of the previous 
subsection, both expressed in the space of spectral parameters. Clearly, the ensemble 
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Figure 5: Constraints from the CMB+LRG7 data set on combinations of the spectral 
parameters (ng ,ag ,r). The black lines denote the constraints computed with the second- 
order slow-roll approximation, starting from flat priors on the HFF parameters at the 
pivot scale (ei, 62, e2e3); the red lines denote the constraints obtained by using directly 
(ns,r,as), but enforcing the second-order consistency conditions for the tensor-to-scalar 
ratio and for the running of the tensor spectral index. 



of models covered by the two parametrisations is the same, and differences might be 
expected only due to different priors. However, the contours of the two runs overlap 
perfectly. This is an indication that in the range of models allowed by the data, the 
expressions for {ns, as, r) in terms of (ei, 62, 6263) are dominated by linear terms, so 
that flat priors on one parameter set are nearly equivalent to flat priors on the other 
one. The run of this section does not bring new information on inflationary models 
but provides a useful self-consistency check. 

Our constraints on the HFFs evaluated at the pivot scale are best summarised 

by 

ei < 0.021 , 62 + 2.7ei = 0.085l°:°i , £263 = 0.0611°;°^^ (95% c.l.). (5.1) 

When €3 is fixed to zero, they reduce to 

ei< 0.011, 62 + 2.761 = 0.040 ±0.026 (95% c.l.). (5.2) 

5.3.3 Flat priors on combinations of the Taylor coefficients of the observ- 
able inflaton potential 

The results of the last two subsections p.3.1|, |5.3.2| should be taken with a grain of salt 
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for two reasons. First, when the running is large, the second slow-roll approximation 
is not necessarily accurate for all models allowed by the data, hence any relation 
between the spectral parameters (InAg, rig, as, r) and the underlying inflationary 
potential ^(0) includes a theoretical uncertainty. Second, for some spectra with a 
large running allowed by the data, there is no guarantee that there exists a single 
underlying inflationary model consistent with these spectra. Since all allowed models 
have r ^ 1, it is clear that they are consistent with the condition ei(A;*) <^ 1 at the 
pivot scale, but on the edges of the observable range, they might be incompatible 
with ei(fc) < 1, i.e., with inflationary expansion. 

In the method consisting in fitting directly the potential V{(f)) with a numerical 
computation of the primordial spectra over the observable range, these caveats are 
avoided by construction. Still, one needs to make an assumption concerning the 
shape of the infiaton potential in this range. In this subsection we assume that 
the observable infiaton potential can approximated by a Taylor expansion of order 
either two or three. Our free parameters are the potential and its derivatives with 
respect to the infiaton field, evaluated when the pivot scale crosses the Hubble 
radius during inflation: Vq = V, Vi = dV/d(f), V2 = (fV/d<j)'^ and = d^V/dcj)^. 
In order to avoid complicated parameter degeneracies, we impose flat priors on the 
combinations (Vl/Vq)^, V2/V0 and V^Vi/Vq-. these parameters are related hnearly to 
the usual 'potential slow-roll parameters' and, although no slow-roll approximation 
is performed here, they remain not too far from linear combinations of (ng, ag, r) 

After running CosmoMC, for each model in our chains we can compute (ng, ag, 
r) defined at the pivot scale directly from the numerical primordial spectra. Hence, 
the results of this run can be readily compared in spectral parameter space with 



those of the previous subsection based on HFF parameters (like in |j78|, ^). This 
comparison is illustrated in Figure |[ As expected, the run with V3 = leads to 
results very similar to those of the power-law model with ag = 0. The reason is 
that a quadratic potential cannot generate large running. This can be proved using 
the second-order slow-roll approximation, which remains accurate for all quadratic 
potentials compatible with the data. Following the same logic, one could expect the 
run with V3 7^ to mimic the results of the running model. This is far from being the 
case, as can be seen in Figure |^ when comparing the black and green contours. The 
model with cubic terms in the potential cannot not reach such large values of |ag| and 
r as the running model. This can be understood in the following way. For models 
with large running and tensors, the Horizon Flow Functions ei(iV^) and e3(A^*) are 
not so small when the pivot scale crosses the Hubble radius, and the flow of equations 
leads ei(iV) to reach order one even within the observable range. So, these models are 
not consistent with any underlying potential described by a third order polynomial 
within the observable range, even without assuming any extrapolation scheme beyond 
this range. The limit between consistent and inconsistent models depends crucially 
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Figure 6: Constraints from the CMB-I-LRG7 data set on combinations of the spectral pa- 
rameters (ns ,as,r) mider various assumptions: (blue) exact computation of the spectrum 
when the potential V{(p) is Taylor-expanded at order two in the observable range; (green) 
same with order three Taylor expansion; (black) spectrum computed with the second-order 
slow-roll approximation, starting from flat priors on the HFFs at the pivot scale (ei, £2, 
£2 £3); (red) same with £3 = prior (and first-order slow-roll approximation). 



on the maximum wavenumber /Cmax at which we require the primordial spectra to 
converge after Hubble crossing (see Appendix [A.3| and the discussion in p4|), in our 
case /cmax = 5Mpc~^. For this value and under the assumption that the observable 
infiaton potential can be described by a third order polynomial, we find that r < 0.15 
(at 95% C.I.), i.e., roughly the same bound as in the power-law model. At the same 
time we get as = —0.32 ± 0.26 (at 95% c.l.): negative running is still favoured, but 
with twice smaller values of |as|- 

Note that the method of this subsection could be iterated at higher-order in 
the Taylor expansion of the observable infiaton. The results of ||7S|] based on older 
data suggest that with a fourth-order expansion, the entire range of (as, r) values 
allowed by the data in the running model could be compatible with the V{(f)) model. 
However, this is at the expense of introducing a significant running of the running /3s 
(this parameter is governed by d'^V / dcj)'^ at the pivot scale). In that case the allowed 



ranges for (ns, as, r) would be even larger than found in section |5.3.1 



The results of the runs based on potential parameters and HFF parameters are 
also interesting to compare in the space of potential parameters {Vi/Vq)'^, V2/V0 and 
VsVi/Vq. In order to perform such a comparison, we re-mapped the HFF param- 
eters into potential slow-roll parameters using the formulae in Appendix |A.2|. The 
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Figure 7: Inflation constraints on combinations of the potential derivatives Vi = dV/dcj), 
V2 = (f'V/d(jP', V3 = d^V/d4>^ under various assumption: (blue) exact computation of 
the spectrum when the potential V{(t)) is Taylor-expanded at order two in the observable 
range; (green) same with order three Taylor expansion; (black) spectrum computed with 
the second-order slow-roll approximation, starting from flat priors on the HFFs at the pivot 
scale (ei, £2, £2 £3), converted here into potential slow-roll parameters using exact formulae; 
(red) same with €3 = prior (and first-order slow-roll approximation). Mp stands for the 
reduced Planck mass. 



results are shown in Figure 0. Again, the power-law model and the second order 
expansion in V^(0) provide similar results, as can be seen from the two-dimensional 
likelihood contours in ((Vi/Vb)^, V2/V0) space. For these models the best constrained 
combination is 
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Ml = -0.023 ± 0.013 , (95% c.l. 



(5.3) 



while 



T^] Mj,< 0.020 



(95% c.l.) 



(5.4) 



When the third derivative of the potential is turned on, the combination VsVi/VJf is 
found to be in the range 



V3V1 
^0 



Ml = O.OUt'^oml , (95% c.l.) 



(5.5) 
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while a reconstruction of the potential from spectral parameters (with running) would 
suggest a larger range 



1/2 



M^p = 0.030l°;°J^ . (95% c.l.) 



(5.6) 



This discrepancy is equivalent to the one for as discussed above. For third-order 
polynomials ^(0), the upper bound on {Vi/VqY is the same as in Eq. ( p.4p but the 
preferred combination of the first two parameter becomes 
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Ml = -0.041 ± 0.021 . 
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Table 3: Mean parameter values and bounds of the central 95%-credible intervals for the 
vanilla+r+as and vanilla+as parametrisations. 



6. Conclusions 

We have provided a detailed update on the present observational status of single- 
field inflation, paying particular attention to the robustness of our results to selection 
of data sets and theoretical uncertainties. WMAP5 data, in combination with the 
measurements of small scale CMB experiments and the halo power spectrum of lumi- 
nous red galaxies from the SDSS data release 7 provide increasingly tight constraints 
on the physics of inflation. After the completion of this work, the WMAP7 analy- 



sis of Komatsu et al ||88[ appeared and reached similar conclusions. Still, some of 
the bounds presented here (based on WMAP5 data) are slightly tighter due to the 
inclusion of small-scale CMB experiments and large scale structure data. 

We find that evidence against the Harrison-Zel'dovich spectrum is mounting, 
with a tilt of the spectrum now favoured at the level of about four standard devia- 
tions. Within the (ns,r)-model which corresponds to most slow- roll inflation models, 
we infer constraints of = 0.962]']oo26 and r < 0.17 (at 95% confidence level), which 
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puts, e.g., chaotic A0^-inflation under severe pressure, but leaves the quadratic po- 
tential model inside the 95% c.l. contours. Generally, the tendency of the data to 
prefer a convex shape of the inflaton potential has increased, though many concave 
models are still viable. We have checked that these conclusions are robust with re- 
spect to cosmological model extensions (such as the inclusion of a non-zero neutrino 
mass) and to theoretical advances in modelling the physics of recombination. 

Tantalisingly, there remain mild indications that the data prefer a negative run- 
ning of the spectral index at roughly two standard deviations - indeed, a fair bit 
of allowed parameter space appears to be in conflict with the predictions of models 
with negligible fourth derivative of the potential in the observable range. We have 
also presented a comparison among different methods in sampling the inflationary 
parameters including running. 

With a multitude of upcoming precise measurements of CMB polarisation and 
smaller scale CMB temperature perturbations by experiments like Planck, ACT, 
SPT, QUIET, SPIDER, PolarBear and EBEX, it is does not seem too optimistic to 
expect that we will soon be able to resolve these issues and take another big step 
towards reconstructing the physics that governed the inflationary era of the Universe. 
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A. Inflation predictions 

We compare infiationary predictions with observations, by adopting a parametrisa- 
tion of the primordial power spectra (PS henceforth) of curvature and tensor pertur- 
bations as: 

.^.,.,,.,(^),!^..(^) 
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where X = S , T stands for scalar and tensor, respectively, fc* is the pivot scale, 
bis = ris - I, biT = nT, 62s = "s, &2T = "t- 



A.l Sampling the spectral parameters 

In Section ^ we have used the parametrisation of Eq. (|A.1|) , sampling directly from the 
logarithm of As (= Vs{K)e^°^), ns and r (= VT{K)e'"'^~'''>^ /VsiK)), with 62X = 0. 
The tensor spectral index has been fixed through the (first-order) consistency condi- 
tion: 

nT = . (A.2) 

In Section |]we have sampled As, ns, r, as- The tensor spectral index has been 
fixed through the (second-order) consistency condition: 

nT = --(2---ns) . (A.3) 

The running of the tensor spectral index has been fixed to: 

T ( T \ 

ot = - - + ns - 1 . (A.4) 



A.2 Sampling the Horizon Flow Functions (HFFs) 

In Section ^ we have derived Vx{k*) ,bix from the Hubble parameter H and the 
horizon flow functions ti (HFF henceforth) evaluated at the pivot scale k^. The 
HFFs are defined as ei = —H/H"^ and Cj+i = ei/{Hei) = (dej/dA^)/ej with i > 1 
and the number of e- folds (dA^ = Hdt) |^^. The analytic slow-roll approximated 
power spectra has been obtained first through the Green's function method (GFM 
henceforth) in Refs. [|^, The coefficients for the scalar spectrum are: 

6so = -2(C + l)ei-Ce2+ (-2C+^-7) e? + - l) el 

+ {-C' - 3C + ^ - 7) e,e2 + {-IC^ + ^ (A.5) 

6si = ns-l = -2ei-e2-2e?-(2C + 3) ei £2 - C ^ (A.6) 

de2 

bs2 = as = -2eie2 - (A.7) 
and for tensors are: 

6to = -2 (C + 1) ei + (-2C + - 7) + (-C^ _ 2C + g - 2) €^€2, (A.8) 

6ti =nT = -2ei - 2el - 2 {C + 1) d €2, (A.9) 

6t2 = Ot = -26162, (A. 10) 

where C = In2 + 7E — 2 ^ —0.7296 (7E is the Euler-Mascheroni constant) and 
d62/dA^ = 6263. A similar structure has been obtained through the method of Com- 
parison Equation 
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For a Klein-Gordon scalar field 0, the potential V and its derivatives are related 
to H and e,- as: 



V 3 



V' ■ (1-f)- 

T/ /\^2 O, £2 I 5£l£2 "^2 ';2£3 

KaiaiMp _ ^ei - Y 3- H g Y2 6~ 



V^^^V^M^ 1 — ^ + ^ 



P _ ^ t + t ^ 2 ^ , £263 2^ 2 

— —3- (^4ei - deie2 + — + 6^62 - eiea 

The use of the HFF parametrisation allows to reconstruct the derivatives of the 
potential without any additional approximation. 

A. 3 Sampling the potential parameters 

In Section ^ we used the CosmoMC inflationary module^ released together with [[7^]. 
The basic principles of this module are the following. 

In CosmoMC, the pivot scale is fixed once and for all, but different parameters 
defining the function V{(f) — (p^,) are passed to CAMB (here 0^, is the value of the 
inflaton field when the pivot scale crosses the Hubble radius; it does not need to be 
formulated explicitly). For each V{(f) — 0*), the module computes the spectra Vs{k), 
V"[{k) within the range [fcmin, /Cmax] = [5 x 10~^,5] Mpc~^ needed by CAMB, imposing 
that aH = when = 0^. So, the code first finds the inflationary attractor solution 
around = 0*, computes if* and normalizes the scale factor so that a* = /c^/if*. 
Then, each mode is integrated numerically for k/aH varying between two adjustable 
ratios: here, 50 and 1/50. The evolution of each scalar/tensor mode is given by 



d?7^ 



+ 



1 d^^s.T 



es,T = (A.12) 



with Tj = J dt/a{t) and zs = a(j)/H for scalars, zt = a for tensors. The code 
integrates this equation starting from the initial condition ^s,t = e~*'^''/A/2A; when 
k/aH = 50, and computes 

Vs = hm , Vt = hm 2"^- • (A.13) 

k<^aH 2%^ ZS, k<^aH nnip zj. 

So, the earliest (latest) time considered in the code is that when kmin/dH = 50 
{kmax/ciH = 1/50), which in the attractor solution uniquely determines extreme 



^http : //wwwlapp . in2p3 . f r/~valkenbu/inf lationH/ 
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values of (0 — 0*) according to some potential. In the module this is translated 
to demanding that aH grows according to the aforementioned ratios: by SO/c^/fcmin 
before = 0*, and by SOfcmax/^* afterwards. Hence, one of the preliminary tasks 
of the module is to find the earliest time. If by then, a unique attractor solution 
for the background field cannot be found within a given accuracy (10% for 0ini), the 
model is rejected. So, the module implicitly assumes that inflation starts at least 
a few e-folds before the present Hubble scale exits the horizon (this assumption is 
relaxed in the second version of the module based on -ff(0) reconstruction ||8^). In 
addition, the module imposes a positive, monotonic potential and an accelerating 
scale factor during the period of interest. As a result of the chosen method, the 
potential is slightly extrapolated beyond the observable window, in order to reach 
the mentioned conditions for the beginning and ending of the numerical integration. 
The range of extrapolation is still very small in comparison with an extrapolation 
over the full duration of inflation after the observable modes have exited the Hubble 
radius. Note that in this approach one does not need to make any assumption about 
reheating and the duration of the radiation era. 
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